Advanced Reactor Modeling with MATLAB by Riccardo Tesser Vincenzo Russo

Author:Riccardo Tesser, Vincenzo Russo
Language: eng
Format: epub
Publisher: De Gruyter
Published: 2020-11-05T12:33:31.671000+00:00

Example 5.1 Tank in series

For the TIS model, implement the following equation varying Θ = 0:10–2:3 and n = 1:100:

Matlab code for the solution of this example and results are reported as follows Figure 5.2:

Figure 5.2: Dimensionless RTD function profile as a function of the dimensionless time and the number of CSTR (left). E(Ï´) versus Ï´ plot for different values of the number of CSTR (right).

% example 3.1

clc, clear t = [0:1e-2:3]' ; n = [1:100]' ; for k = 1:100 E(:,k) = (k.*(k.*t).^(k-1)).*exp(-k.*t)./factorial((k-1)) ; end figure(1) subplot(1,2,1) surf(n,t,E,'EdgeColor','none') xlabel('\it heta
m[-]'), ylabel('\it n
m[-]'), zlabel('\it E_{ heta}
m[-]') subplot(1,2,2) plot(t,E(:,2),'-r',t,E(:,4),'-k',t,E(:,10),'-b',t,E(:,100),'-g') xlabel('\it heta
m[-]'), ylabel('\it E_{ heta}
m[-]') legend('\it n
m = 2','\it n
m = 4','\it n
m = 10','\it n
m = 100') grid on

Higher number of CSTR leads to a sharper peak. With an infinite number, in theory the output would give a δ-Dirac function, characteristic of a plug flow.

Download

Copyright Disclaimer:
This site does not store any files on its server. We only index and link to content provided by other sites. Please contact the content providers to delete copyright contents if any and email us, we'll remove relevant links or contents immediately.